Relaxation dynamics of an elastic string in random media 
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We investigate numerically the relaxation dynamics of an elastic string in two-dimensional random 
media by thermal fluctuations starting from a flat configuration. Measuring spatial fiuctuations of 
its mean position, we find that the correlation length grows in time asymptotically as ^ ~ (Ini)^''*. 
This implies that the relaxation dynamics is driven by thermal activations over random energy 
barriers which scale as -E's(^) ~ with a length scale £. Numerical data strongly suggest that the 
energy barrier exponent x is identical to the energy fluctuation exponent x = 1/3- We also find 
that there exists a long transient regime, where the correlation length follows a power-law dynamics 
as 5 ~ t^^^ with a nonuniversal dynamic exponent z. The origin of the transient scaling behavior is 
discussed in the context of the relaxation dynamics on finite ramified clusters of disorder. 
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Interaction and quenched disorder are essential ingre- 
dients in condensed matter physics. Many-body systems 
may undergo a phase transition as a cooperative phe- 
nomenon mediated by interaction. When quenched dis- 
order comes into play, there may emerge a glass phase in 
which degrees of freedom are pinned by random impuri- 
ties and a slow relaxation dynamics appears. An elastic 
string in random media is one of the simplest systems 
where the interplay between interaction and quenched 
disorder yields a nontrivial effect [l|, Q ■ This has been 
studied extensively in literatures since it is relevant to 
many interesting physical systems such as a growing in- 
terface 3 , a domain wall in random magnets 0, Q , and 
a magnetic flux line in superconductors [6| . 

In low temperatures, one may approximate an elastic 
string as an elastically coupled directed polymer where 
no overhang is allowed. Then it can be described by a 
single valued function (m), where x £ TV^ and u £ TZ are 
the transverse and the longitudinal coordinates to the 
polymer direction, respectively, in a (d -f 1) dimensional 
space. The energy of a polymer of length L in a configu- 
ration x{u) with < u < L is given by the Hamiltonian 
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The first term accounts for an elastic tension and the 
second term V{u, x) is a random pinning potential with 
short-range correlations. 

Equilibrium properties of the directed polymer in ran- 
dom media (DPRM) are rather well understood. The 
tension favors a fiat state, while thermal fiuctuations and 
the disorder potential favor a rough state. The compe- 
tition between them leads to the scaling law |Ax| ^ L'' 
for the transverse fiuctuation (interface roughness) and 
AE ~ for the (free) energy fluctuation. The quenched 
disorder is relevant for d < 2, and the polymer is in a 



super-rough phase (C > 1/2) at all temperatures. Espe- 
cially for d = I, the scaling exponents are known exactly 
as (id = 2/3 and xiD = 1/3 4]. For d > 2, it is believed 
that there is a transition from a super-rough phase into 
a normal-rough phase {( = 1/2) as the temperature T 
increases In the latter, the thermal fluctuations dom- 
inate over the disorder fiuctuations, while vice versa in 
the former. 

When a polymer is in a nonequilibrium state, e.g., a 
fiat configuration, it will relax to the equilibrium rough 
state. Without disorder or in the normal-rough phase 
with disorder for d > 2, the elastic polymer equilibrates 
diffusively. The correlation length ^ in the longitudinal 
u direction grows algebraically in time as f ~ with 
the dynamic exponent Zq — 2. In the presence of the 
quenched disorder, one expects a slower relaxation be- 
cause random impurities tend to trap the polymer into 
metastable states in local energy valleys. Upon equili- 
bration, the polymer has to overcome energy barriers Eb 
separating those valleys through thermal fluctuations to 
approach the true equilibrium state. It is believed that 
the energy barrier height scales as Eb {t} ^ \n a region 
with linear size £. 

Thermal activations allow the correlated polymer seg- 
ment of length ^ to overcome the energy barriers in a time 
scale ~ QEB{i)/T ^ Then, it follows that the correlation 
length ^ grows as 

m - (rini)i/x, (2) 

with the universal energy barrier exponent x, indepen- 
dent of the disorder strength and the temperature 
Assuming that there is only a single relevant energy scale 
in this system, the exponent yshould be equal to the en- 
ergy fluctuation exponent x \M- This conjecture is sup- 
ported at least in low dimensions • 

However, even for d = 1, there is a long-standing con- 
troversy on the scaling law of Eq. ([2]). Numerical simu- 
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lation study [8| reports a signature of the expected log- 
arithmic scaling, but only after a long and clean inter- 
mediate power-law scaling regime where ^ ~ t^/^ with 
a nonuniversal dynamic exponent z, whose origin is not 
clear. Moreover, the scaling exponent associated with 
the logarithmic scaling seems different from the conjec- 
tured value of X ^ X ~ 1/3 @]- There is also a recent 
claim of X = d/2 based on the droplet theory [§]. Besides, 
there are many numerical works in the context of domain 
wall coarsening dynamics in two dimensional random fer- 
romagnets |5i], which seem to support the nonuniversal 
power-law scaling without any signature of the asymp- 
totic logarithmic scaling 13, U | ■ These intriguing results 
are also left unexplained. 

In this work, we study the relaxation dynamics of the 
DPRM in (1 -I- 1) dimensional lattices with extensive nu- 
merical simulations and a scaling theory approach. Mea- 
sured are the spatial fluctuations of the mean position of 
the polymer, from which we derive the dynamic scaling 
behavior of the correlation length £^{t) through a simple 
scaling hypothesis. The purpose of this study is to settle 
down the controversy by providing a decisive numerical 
evidence on the asymptotic relaxation dynamics of the 
DPRM for d = 1. In addition, we suggest a reasonable 
scenario for the origin of the transient power-law scaling 
regime. 

We consider a discrete model for the DPRM in the 
45-degree- rotated square lattices of size L x M. Each 
lattice site is represented as {i,x) with the longitudinal 
coordinate i — 0, ■ ■ ■ , L — I and the transverse coordi- 
nate X = —M/2 -t- 1, • • • , M/2 with the constraint i = x 
in modulo 2. Assigned to bonds are quenched disorder 
variables J which are distributed independently and ran- 
domly according to a probability density function p{J). 
The polymer of length L is placed along the bonds and 
directed in the longitudinal direction without any back 
bending. Then its configuration is described by the fluc- 
tuating variables {x{i)} with the solid-on-solid (SOS) 
constraint of \x{i) — x{i ± 1)| = 1 for all i. We adopt 
the periodic boundary condition in the longitudinal di- 
rection, i.e., x{L) — x{0). The transverse size M is taken 
to be large enough (A/ = 4096 ~ 8192) to avoid any 
possible interference. 

The polymer energy is given by the lattice Hamiltonian 



W = ^ J(i, xii);i + 1, x{i + 1)) , 
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where J{i,x;i + l,x') denotes the disorder strength of 
the bond between neighboring sites (i, x) and {i + l,x'). 
In this study, we consider the uniform distribution in 
the range — 1 < J < 1 and the bimodal distribution 
p{J) = fS{J + 1/2) + (1 - f)S{J - 1/2) with a model 
parameter /. It turns out that both cases lead to the 
same conclusion. 

We start with the flat configuration with {x{i) — i 



mod 2} at t = and study its relaxation dynamics to- 
ward the equilibrium state at temperature T. We adopt 
the Glauber Monte Carlo dynamics: First select a site i 
at random and try to flip x(i) x(i) ± 2 with prob- 
ability 1/2, respectively. Unless the trial violates the 
SOS constraint, it is accepted with the probabihty of 
min[l, e~^^/-^] where AE is the energy change. The 
time is incremented by one unit after L such trials. 

We focus on the spatial dispersion of the mean position 
X = '^iX{i)/L as function of elapsed time t, which is 
given as 
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Here (•)t and [•]/) denote the thermal and disorder aver- 
age, respectively. Let £^{t) be the characteristic correla- 
tion length of the polymer in the longitudinal direction. 
At time t, each segment of size ^{t) equilibrates with 
transverse displacement of the order 6x ~ with the 
roughness exponent C. When ^ <IC L, each segment is 
independent and the total displacement is given by 



(Ax) 
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Utilizing this relation, we can derive the correlation 
length ^(t) from the ensemble-averaged global quantity 
{Ax)'^{t), which usually bears a better statistics than the 
distance-dependent correlation function of the transverse 
displacement. 

When the polymer fully equilibrates, i.e., ^(t = r) ~ L 
with the relaxation time r = t(L), the polymer as a 
whole (the mean position) starts to diffuse normally in 
the transverse direction. One expects that (Ax)^ grows 
linearly in time scaled by r as 



(Ax)^ 



(6) 



Without disorder, the motion of the polymer is governed 
by the linear Edward- Wilkinson (EW) equation [I4I . The 
EW class is characterized by = 1/2 and f ~ t^/^, i.e., 
T ~ L^. In this case, the dynamic behaviors before and 
after equilibration, Eqs. (O and ([5]), follow the same scal- 
ing law (Ax)^ ~ t/L at all t. Indeed, this coincides with 
the exact solution of the EW equation. However, in gen- 
eral with disorder, these two scaling laws are distinct. 

We have performed extensive Monte Carlo simulations 
to examine the scaling property of (Ax)^. Figure [1] shows 
a plot of the numerical data with the bimodal disorder 
distribution with / = 0.25 at T = 1.0. As the scaling 
form predicts in Eq. ([S]), (Ax)^ is inversely proportional 
to L for t < t{L), so it is convenient to plot L(Ax)^ 
versus t, where all curves with different L collapse into 
one scaling curve for t < t{L) and then start to deviate 
and show the finite size effects given by Eq. 

We find that there exist four distinct regimes: (I) 
For t < to(~ 10^), the polymer moves diffusively as 
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FIG. 1: Monte Carlo simulation data for the system with the 
bimodal disorder distribution with / — 0.25 and at T = 1.0. 
Different symbols represent data for different polymer lengths 
L. The data are averaged over Ns ~ 5000 disorder samples. 



L{Ax)^ ^ t. In this regime, the polymer behaves as 
in the EW class since it does not feel the disorder pin- 
ning as yet. (II) For to < t < tc(~ 10'*~^), the poly- 
mer is affected by the disorder and exhibits a power-law 
scaling behavior as L(Aa;)^ ~ with a nonuniversal ex- 
ponent a. The crossover time tc is very large but finite 
and independent of L, which implies that this power-law 
scaling is transient. (Ill) For tc < t < t{L), there is a 
continuous downward curvature in the plot suggesting a 
possible logarithmic scaling with L(Aa;)^ ~ (Ini)"^. The 
crossover time t(L) increases indefinitely with L, which 
implies that this regime should be the true asymptotic 
scaling regime. (IV) For t > t{L), the polymer displays 
a diffusive motion with a size-dependent diffusion ampli- 
tude as in Eq. ([6]). 

We investigate the scaling behavior in each regime 
quantitatively. Useful are the effective exponents de- 
fined as aeff{t) = dln(L(Aa;)^)/dlnt and (j)eff{t) = 
dln(L(Aa;)2)/dlnln<. If L{Axf - as in the regime 




FIG. 2: Effective exponent plots for the data given in Fig. \T\ 
Our estimates are a = 0.63±0.03 in regime II and (ji = 6.8±0.5 
in regime III. 
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FIG. 3: Effective exponent plots with the uniform disorder 
distributions with the temperature T = 0.25 in (a) and T = 
0.75 in (b), and the bimodal disorder distribution with / = 0.1 
and r = 1 in (c) and / = 0.5 and T = 1 in (d). The symbols 
have the same attribute as in Fig. [T] The dashed lines are 
drawn at cj)eff = 7. 



II, one would obtain that aef f{t) = a and (j)ef / (t) = 
alnt. On the other hand, if L(Ax)^ ~ (Int)'^ as in the 
regime III, one would obtain that aeff(t) = (j)/\nt and 

<l>eff{t) = 0. 

We plot the effective exponents in Figs. [2] (a) and (b). 
In (a), one can clearly see a plateau at aeff — 0.63 ± 
0.03 for 10^ < t < 10'*(= tc) (regime II). For t > tc, 
it continuously decreases in the regime III before hiking 
up in the regime IV as expected. In (b), there appears 
a plateau at 0e// = 6.8 ± 0.5 in the regime III, which 
widens as L increases. These numerical evidences lead to 
a definitive conclusion that the asymptotic motion of the 
polymer follows the logarithmic (not power-law) scaling 
as 



L{Axf{t) - (Int)" 



(7) 



with the exponent (f> — 6.8 ± 0.5. In terms of the corre- 
lation length, we find, using Eq. 



1/x 



(8) 



where 1/x = </'/(l + 2C) = 2.9 ±0.2. The estimated value 
of X = 0.34 ± 0.03 clearly favors the conjecture value of 
X = X = 1/3 and invahdates the recent claim of x = 1/2 
(or equivalently (j) = 14/3). 

In order to examine the universality of the logarith- 
mic scaling behavior, we have performed the simulations 
with different values of / in the case of the bimodal dis- 
order distribution and also in the case of the uniform 
disorder distribution. The effective exponents are then 
presented in Fig. [3l The plateaus around 0e// = 7.0 
shown in Figs.[3l^a), (c), and (d) confirm the universality 
of the logarithmic scaling as well as the universality of 
its exponent x- The plot in Fig. ^h) seems to be in- 
compatible with (j) = 7. Comparing Figs. [3Ka) and (b). 
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one can notice that the finite-size behavior (regime IV) 
sets in earher at the higher temperature. So the loga- 
rithmic scahng regime is observed only when L > 256 
in (b), while it is already evident at L = 64 in the other 
cases. This suggests that one would need larger polymers 
at higher temperatures. 

It is puzzling why there exists the extremely long tran- 
sient regime II where the polymer relaxation seems to fol- 
low a power-law scaling such as L{Ax)^ ^ or ^ ~ t^^^ 
with z — {l + 2(^)/a. Moreover the exponent a is nonuni- 
versal and varies with the disorder strength and the tem- 
perature. Such a transient behavior was also reported in 
Ref. [1, [3], but its origin has never been explored. We 
suggest one reasonable scenario as below. 

It is convenient to consider the bimodal disorder dis- 
tribution. When / < 1/2, the energetically favorable 
bonds with J — —1 may play the role of local pinning 
centers for the polymer [ij]. As the polymer consid- 
ered here is directed, relevant are the directed percolation 
clusters of the pinning bonds. These clusters are rami- 
fied but finite in size, as the bond density / is smaller 
than the directed per colation threshold fc ~ 0.6449 in 
the square lattice [15|. The characteristic size and the 
mean distance between them are denoted by and ^i, 
respectively. After the initial diffusive motion, polymer 
segments are trapped by those clusters independently as 
long as the correlation length is smaller than the cluster 
size < lo)- 

The pinning mechanism in fractal-like ramified lattices 
is different from that in bulks. The energy barrier height 
in such lattices is shown to scale logarithmically with a 
length scale £ as Eb {() — Eolni with an universal con- 
stant Eq depending only on the ramification degree [l3| ■ 
Then, the time scale associated with the thermal acti- 
vation of the correlated segment of length ^ is given by 

power-law growth 

of the correlation length as ^ ^ ^^^^ with the 

nonuniversal dynamic exponent z = 2 + Eo/T. The tem- 
perature dependence seems consistent with our numerical 
estimates for z (not shown here). When ^ exceeds Iq, the 
polymer segments are pinned by a few pinning clusters. If 
£, > h, then the polymer starts to be pinned collectively, 
and the transient power-law scaling behavior crosses over 
to the asymptotic logarithmic scaling behavior. 

Finally, we add one remark on the domain coarsening 
dynamics in the two-dimensional random ferromagnets. 
When the system is quenched well below an ordering tem- 
perature from a disordered state, the characteristic size R 
of ordered domains increases and the domain wall motion 
may be described by the DPRM. Hence, it is natural to 
expect that R{t) - ^^-C _ (lni)(2-C)/x Surprisingly, 
recent high accuracy numerical simulation studies re por t 
that R{t) ^ t^/^ with a nonuniversal exponent z 
Our result suggests that those behaviors may be due to 
the pinning of domain walls by finite pinning clusters in 
the transient regime. 



In summary, we have investigated numerically the re- 
laxation dynamics of the DPRM. The numerical data 
show unambiguously that the correlation length grows 
as ^ ~ t^^^ in the transient regime and then ^ ~ (Int)^/^ 
in the asymptotic regime. The transient behavior is orig- 
inated from the pinning independently by local ramified 
impurity clusters. The asymptotic logarithmic scaling 
is compatible with the scaling picture that the energy 
barrier height scales in the same way as the energy fluc- 
tuations with X = X = 1/3. Implication on the domain 
coarsening dynamics is also discussed. 
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